Critical interfaces in the random-bond Potts model 
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We study geometrical properties of interfaces in tlie random-temperature g-states Potts model as 
an example of a conformal field theory weakly perturbed by quenched disorder. Using conformal 
perturbation theory in g — 2 we compute the fractal dimension of Fortuin Kasteleyn domain walls. 
We also compute it numerically both via the Wolff cluster algorithm for g = 3 and via transfer-matrix 
evaluations. We also obtain numerical results for the fractal dimension of spin clusters interfaces 
for 5 = 3. These are found numerically consistent with the duality k'^'^^'^kF^ = 16 as expressed in 
putative SLE parameters. 



The discovery of the Schramm Lowner Evolution 
(SLE) has strongly revived interest in geometrical prop- 
erties of interfaces in 2-dimensional statistical physics. 
SLE provides a rigorous classification, with a single pa- 
rameter K, of probability measures on non-crossing ran- 
dom fractal curves, which satisfy both conformal invari- 
ance and the domain Markov property [T] . Interfaces and 
similar geometric objects defined in pure 2-dimensional 
critical models are conjectured, and in some cases proven, 
to satisfy both requirements in the continuum limit. 
SLE hence describes such diverse systems as percolation 
K — 6, self-avoiding walks k — 4/3, loop-erased random 
walks K = 2 and level lines of height models k = 4 [2j. 
It applies to the Ising and 3-states Potts interfaces, both 
for spin clusters {k — 3 and k = 10/3, respectively) and 
the dual Fortuin Kasteleyn (FK) clusters {k — 16/3 and 
K = 24/5), with a duality k ^ k' = 16/k. While these 
models have been described, prior to SLE, using con- 
formal field theory (CFT), SLE bridges the gap between 
the algebraic approach of CFT and the geometry of inter- 
faces. The present connections between SLE and CFT [2] 
focus on boundary-condition changing operators, which 
generate the curves. They give df — 1 + k/8 for the frac- 
tal dimension of the interface, i.e. the hull of the SLE 
trace. Extensions beyond non-minimal CFT [3] are rare. 

A tantalizing question is whether CFT methods and 
SLE help to understand a broader class of scale-invariant 
2d complex systems, such as systems with quenched dis- 
order or far from equilibrium. Numerical studies indicate 
that zero-vorticity lines in 2D-turbulence [4] and domain 
walls in spin-glasses [5] may be described by SLE. These 
examples are "far" from any pure CFT, thus the situa- 
tion may be more favorable for models which are "weak 
perturbations" of a known CFT. This is e.g. the case for 
the g-states Potts model, perturbed by quenched random 
bond (i.e., temperature) disorder, known to exhibit a sta- 
ble weak-disorder fixed point for q > 2, perturbatively ac- 
cessible in a g — 2 expansion. This has been studied using 
conformal perturbation theory and transfer-matrix 
methods 8 . However, geometric properties of interfaces 
which are crucial for future comparison to SLE were to 
our knowledge not investigated. 



The aim of this Letter is to present results for the frac- 
tal dimension of domain walls in the random-temperature 
Potts model. These are obtained by analytical calcu- 
lation using conformal perturbation theory inspired by 
[6l[7], and from two types of large-scale numerics: Monte 
Carlo simulations using the efficient Wolf-algorithm [TD] , 
which allow to keep track of both spin and FK clusters 
in the same simulation, and transfer-matrix calculations, 
whose advantage is to make close contact with CFT. The 
results of all three methods agree nicely. 

Let us recall the definition of the model: In terms of 
the spin variables cr^ — {I,-- - ,q} at lattice site i, the 
partition function of the g-states Potts model is 

where the sum runs over nearest-neighbor bonds (i, j). 
The last expression is the spin-cluster expansion, not- 
ing 1 — Pij = gy expanding in Pij, it can be 
rewritten in terms of the Fortuin-Kasteleyn (FK) clus- 
ters, composed by placing a bond between neighboring 
sites with probability pij . The pure ferromagnetic model 
has Jij = J > 0, Pij = p, while in the disordered one the 
Jij are chosen as i.i.d. random variables. The partition 
function in the FK representation is (up to a prefactor) 

Z^^piei(l-p)l^lglieil , (1) 

e 

for the pure model, with a straightforward generalization 
to the random case. Here Q runs over all clusters (i.e. 
domains connected by the above placed bonds), \G\ is the 
number of bonds, \Q\ the not placed bonds, and \ \Q\ \ the 
number of connected components. The partition sum Q 
allows to define the Potts model with non-integer q> 0. 
For the pure model it has a continuous phase transition 
for < g < 4, which becomes first order for q > 4. 

Our analytical calculation focuses on weak disorder, 
where the Jij = J + 5 Jij are Gaussian random vari- 
ables of variance f3'^6Jfj — go and <C /3J. Near the 
critical temperature of the pure model, the continuum 
limit of the random Potts model can be written j6j |7] as 
'H = '^purc + Jg£{z)SJ{z) where Jg= f d^z and /3Hpurc is 
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the action of the pure q-state Potts model, which at crit- 
icality can be identified with its conformal field theory, 
or the 0{N = y/q) model [12]. We use the Coulomb-gas 
representation of the latter. The coupling constant p is 
related to q via y/g = 2cos(7r/(2p)), so that p — 2 is 
Ising and p = 3 is 3-state Potts. The second term in Ti, 
is the deviation from the pure critical point induced by 
the disorder, where e{z) is the energy density operator 
of the pure model. To average over disorder, the n-times 
replicated action is taken: 

n n 
1 J z 



a=l 
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where everywhere below we use the shorthand notation 
^{z) to denote ^{z), where z = {z,z). The diago- 
nal term e°'{z)'^ is perturbatively less relevant [6] than 
£"(z)£^(z), whose dimension is 4:A^, i.e., four times the 
dimension of the holomorphic part of the energy primary 
field e{z) = <I>i2(z), with Ag = 2(2^-1) ' -^^"^ ^^^^^ 
model, p — q — 2, and disorder is marginally irrelevant, 
whereas for the 3-states Potts model p = g = 3 it is 
relevant. Since the Coulomb gas is defined for all p, we 
can perturbatively expand around the Ising model [S]. 
This expansion is conceptually the same as for the 
model, except that Feynman diagrams are evaluated us- 
ing the unperturbed CFT (with averages denoted (...)o). 
We keep the perturbed system on its critical manifold, 
s.t. only the renormalization of the disorder go is left to 
consider, with a correction to second order O(5o)- 

= 4(n - 2) Efc^, sHz)e'^{z) {s{z)e{z')), + ■■■ . (2) 

Using the unperturbed average (£(z)£(z'))g = |^_ J|4a, 
one obtains the renormalized disorder gL^^^~'^ = 5o + 
47r(n — 2)ff§ 2-4A^ ' ^ being the infrared cutoff, and the 
/3-function (for q> 2) [B]: 

LaL5= (2-4A,)g + 4^(n-2)52 + ... (3) 

At n = 0, (3{g) has an infrared fixed point at g* — ^^l^' 
which determines the low-energy behavior of the random 
model. Conformal symmetry is expected to be restored 
at g* . To date, this method has been employed to cal- 
culate the scaling dimension of the energy density £ and 
of the spin cr, to two- and three-loop order in [6l [7] and 
[7] respectively. The multi-scaling properties of spin-spin 
correlation function has been determined in [3]. 

Here we focus on geometrical properties, hence on the 
operator '^io{z) which measures [O] the passage of one 
critical curve at point z. Indeed, for the pure model. 



gives 



the correlation function ($io(2;)$io(0))o = \z\~'^^^'' 
the probability that two points lie at the perimeter of 
the same FK cluster, from which one obtains the fractal 
dimension of FK domain walls d^^'^"™ = 2 — 2Aio, i.e. 



df^ = 8/5 for g = 3. Here we compute the corresponding 
probability for the disordered system. A crucial question 
is whether $10(2;) is still the "curve-detecting" operator 
in the disordered system. This is true at the "critical 
dimension" p ^ q = 2. Increasing p deforms the operator 
adiabatically. Since the latter is a physical observable, it 
is an eigenoperator of the RG. We must check if there is 
an operator at p — 2 which (i) has the same dimension 
as <&ijOi E^nd (ii) appears in the sub-algebra generated by 
<i>i^o ^iid $1^2- If such an operator exists, it mixes with 
<i>i^Oi and the curve-detecting operator will be one of the 
eigenoperators of the RG flow involving <I>io. We checked 
the absence of such an operator: thus, at least for small 
p — 2, <i>io is the curve-detecting operator. 

We now sketch the calculation of the scaling dimension 
of $10, for details see [IS]. There is no contribution to or- 
der go, since contracting the disorder operator J^b^c ^^^"^ 
with <i>io in, say, replica a, leaves one e in replica b ^ a, 
thus is not proportional to <i>°Q. At second order, con- 
tracting two disorder vertices with <i>"g(zi) gives 



*?o(^i) 



E 

b^c 



s''{z,)s%Z3) 



and projecting onto <i>°Q(zi). Contracting using 
(£'^(z2)£'^(2;3))g — 5^'^\z2 — zal^"''^!^ to eliminate rephcas 
not equal to a we obtain ^1o{zi)e'^^'^ {z2)£'^^°' {z^) , which, 
projected onto '^1q{zi) yields 
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{£{z2)£{zz))a 



2, ^v^ -^,,3— --^-'^ (4) 
X {^iQ{z^)e{z2)e{z:i)\^w{zi)) , 



where the OPE coefficient {f^iQ{zi)e{z2)£{z^)\'^iQ{zi)) := 
limR (i'io(2i)g(22)e(23)i'io(J^))o rp, . , intcCTal 

limfl^OO (*io(2i)$io(i?^))o ■ lULLgldl 

^ arc computed using Coulomb gas techniques )llj . 
One 2d integration, over one angle and one scale, is easy, 
and gives a pole in l/{p — 2). One 2d integral over say Z2 
is left, but we also need a screening charge V+ to get the 
4-point function in ^ . We evaluated this integral in the 
marginal dimension, i.e. for p = 2 (Ising) by analytical 
techniques 7 , and numerically [l6 . The result is 

($10(21) £(2:2) £(^3) *io(zi)) {e{z2)e{z3)) 



= -7.0710L'^"^'^i^(l - 2A12) 



-1 



(5) 



Inserting the fixed-point value g* from above gives 

dimi($io) = -2Aio + ii^^ X 7.071 -|-H0.01433. 
This leads to the fractal dimension of FK domain walls: 

dimL($io) = 1.61433 . (6) 

Let us note a few additional peculiar features which come 
out of the calculation [Tn]. The 4-point function 

G(u):= lim |z4|^'^"('i>io(0)£(l)£(u)a>io(^4)) (7) 
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FIG. 1: Fractal dimension of FK and spin clusters, both for 
the pure (green) and disordered (red) system, using the Wolff 
algorithm. 



at p = 2, i.e., for the Ising model is 
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G is the Meijer G-function, which has a logarithmic di- 
vergence at M = 0, 
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dropping regular terms. The structure of the result and 
the logarithmic divergence remain valid for larger values 
of p, with the parameters replaced by rational functions 
of p. This behavior is consistent with the appearance of 
operators of canonical dimensions 1 and (different from 
the identity) in the OPE of e with <I>i g as discussed in a 
similar case in |17j . Logarithms are known to appear for 
operators on the boundary of the Kac table [18 and in 
disordered systems 

Even more surprisingly, we attempted to perform the 
same calculation for the fractal dimension of spin inter- 
faces, using the operator $01 as curve detector. There 
the equivalent of ([t]) does not exist: At least we were not 
able to construct a 4-point function, which satisfies the 
differential equation induced by the 0-vector condition 
associated to e = <I>i2 at level 2, which is unique- valued, 
and reproduces the correct OPE in the limit of w — *■ 1. 

We now discuss our numerical results. For the Monte 
Carlo simulations we use the Wolff cluster algorithm 
[To] . It consists in randomly choosing a spin, and then 
joining with probability p nearest-neighbor spins to be- 
long to the same cluster. This procedure is repeated 
until no nearest neighbor can be joined anymore, p 
itself is a quenched random variable on the edges of 



the lattice, taken from the symmetric bimodal distribu- 
tion {pi, P2} = {1 - exp(-/3cJi), 1 - exp(-/3cJ2)}- The 
choice P1P2 = g ensures that the system is at its critical 
point [H]. We use Ji = J2 for the pure, and J1/J2 — 10 
for the random-bond disordered system. Once the clus- 
ter is constructed, we assign to it with equal probability 
one of the q colors, as long as this is consistent with the 
boundary conditions (BC) discussed below. Wolff has 
shown [in] that this algorithm produces the correct sta- 
tistical weight of the Potts model, is ergodic, and that 
the critical slowing down is reduced compared to ordi- 
nary Monte Carlo algorithms. In addition, it allows to 
track both spin and FK clusters. In the simulations, we 
imposed boundary conditions which create a domain wall 
which span the lattice as was already done by Gamsa 
and Cardy [T5]. We considered various types of con- 
formally invariant boundary conditions: "fluctuating" 
(a/a), "fixed" (a/b), and "free" (a/free), all of them giv- 
ing the same result in the large-size limit. 

We measured the fractal dimension from the average 
length I of the domain wall as a function of the linear size 
L of the lattice (l) ~ L'^'', where (• • • ) denotes the ther- 
mal average and the disorder average [2S] . The results 
presented here are obtained with a thermal average over 
~ 10*^ X T for the pure system and a disorder average 
over ~ 10^ configurations for the disordered system, r 
is the autocorrelation time which was first determined 
for each size, see [TB] for details. Our simulations show 
that for the pure system, all these domain walls have 
asymptotically the same fractal dimension, with the ex- 
ception of the common domain wall for fixed BC, which 
has dimension one. In Fig. [T] we plot the effective fractal 
dimension versus L. As L — > 00 the fractal dimensions 
of the pure system converge to the values predicted by 
conformal field theory, p'" = y| , and df^ = | , corrob- 
orating partial results by Gamsa and Cardy |T5| . Our 
estimate from all BC, extrapolated to an infinite system 
gives dfP'" = 1.416 ±0.002, and dj^ = 1.599 ±0.002. For 
the disordered system we find 



jspin 



1.401 ±0.003 



1.614 ±0.003 (8) 



This is in excellent agreement with our analytical re- 
sult ([6]). We have also checked the SLE duality relation 
kk' = 16. Using = 1 + f , we find for the pure sys- 
tem k"P'" = 3.328 ± 0.016, k^k = 4.792 ± 0.016, and 
^spin^FK _ 15.95 ± 0.13. For the disordered system, we 



find k"P'" 



K 



FK 



4.912 ± 0.024, and 



^spin^FK 



-- 3.208 ± 0.024, 
15.76 ±0.20. 

In the transfer-matrix approach, we studied the FK 
clusters in the equivalent loop formulation [12]. The 
loops are defined on the medial lattice as the external and 
internal hulls of the FK clusters. The random bonds were 
again drawn from a bimodal distribution, with an equal 
number of strong and weak bonds [H]. The strength of 
the disorder is conveniently characterized by the param- 
eter s PO] defined by J1/J2 = ln(l + s^g)/ ln(l + ^g/s). 
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FIG. 2: Effective central cfiarges c(L) and fractal dimensions 
df^(L) versus disorder strenght s, for s = 3.5,4.0,4.5. The 
linear interpolation is only a guide to the eye. Each curve has 
been normalized so that finite-size effects are absent for s = 1 
(no disorder). Horizontal lines give corresponding exact (resp. 
perturbative) values for the pure (resp. disordered) system. 



where the quantities on the r.h.s. are self-averaging. Only 
the first two cumulants contribute significantly. CFT pre- 
dicts [23] that In(Gj)" = ^ + il H , where the 

n-dependent conformal weights Xj are related to the de- 
sired (multi)fractal dimensions via dj = 2— Xj. For n = 1, 
we have d2 = df^ defined above; c?4 gives the dimension 
of "red bonds" (whose removal disconnects a cluster) . As 
seen from the right panel of Fig. |2] the effective values of 
dj depend strongly on s, so accuracy for s* is important 
PP] . Using s* = 4.0±0.3 we estimate ^2 = 1.615±0.002, 
in excellent agreement with ^ and ([s]). 

To conclude, our analytical and numerical results for 
the fractal dimension of FK domain walls agree well. 
Fractal dimensions of spin interfaces have been deter- 
mined from numerics and seem in agreement with the 
duality relation suggested by SLE. Pending questions un- 
der investigation are possible multiscaling, the fractal di- 
mensions of spin interfaces and SLE type observables. 

We thank M. Bauer, D. Bernard, VI. Dotsenko, A. 
Ludwig, T. Quella and P. Wiegmann for valuable discus- 
sions. Supported by ANR (05-BLAN-0099-01). 



with q = 3. For a given fixed realization of the random 
bonds on long cylinders of length M = 10^ and circum- 
ference L — 4, 6,... ,12 lattice spacings (for the medial 
lattice) the corresponding free energy fj{L), normalized 
per lattice site, was computed from the leading Lyapunov 
exponent of the corresponding product of random trans- 
fer matrices. The transfer direction was taken axial with 
respect to the medial lattice (hence diagonal with respect 
to the original square lattice supporting the Potts spins) 
[TO] . Three different topological sectors were considered, 
corresponding to enforcing j = 0,2,4 loop segments to 
propagate along the length direction of the cylinder. The 
fluctuations of these free energies were studied by aver- 
aging over at least M' = 10^ independent cylinders. 

Conformal field theory predicts [22] that fo{L) — 
/o(oo) — + pr + ■ ■ • , where c is the effective central 
charge and A a non-universal constant. Applying this 
to three consecutive L gives estimates c{L — 4, L — 2, i) 
shown in the left panel of Fig. |2] The fixed-point value 
s* of the disorder strength corresponds to the locus of 
the maximum of c, and is estimated as s* = 4.0 ± 0.3 
(using also data not shown here), improving on the 
value s* — 3.5 ± 0.5 reported earlier [20]. The effec- 
tive central charge of the disordered model is estimated 
as c(s*) = 0.8024 ± 0.0003, in excellent agreement with 
the three-loop perturbative result [7j c ~ 0.8025. 

Correlation functions Gj{M) are defined as the prob- 
ability of having j loop segments propagate over a dis- 
tance M along the cylinder axis without joining up. They 
are related to the free energy gaps through Afj{L) = 
fj{L) — fo{L) = -^\nGj{M). Their disorder-averaged 
n'th moment can be extracted from the cumulant expan- 
sion gj In(Gj)" = nEG7 + \n^{\TiGj -WOjY + • • • , 
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